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Using the formalism of the Khalatnikov potential, we derive exact general formulae for the en- 
QO ■ tropy flow dS/dy, where y is the rapidity, as a function of temperature for the (1+1) relativistic 

hydrodynamics of a perfect fluid. We study in particular flows dominated by a sufficiently long 
hydrodynamic evolution, and provide an explicit analytical solution for dS/dy. We discuss the the- 
oretical implications of our general formulae and some phenomenological applications for heavy-ion 
^ ■ collisions. 
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I. INTRODUCTION 



There is an accumulating evidence that hydrodynamics may be relevant for the description of the medium created 
in high-energy heavy ion collisions [1]. Indeed, experimental measurements such as the elliptic flow [2] shows the 
existence of a collective effect on the produced particles which can be described in terms of a motion of the fluid. 
1 " 1 ' More precisely, numerical simulations of the hydrodynamic equations describe quite well the distribution of low-pj_ 
particles [1], with an equation of state close to that of a "perfect fluid" with a rather low viscosity. On the other 
■ hand, it seems useful to discuss a simplified picture [3, 4], which can be qualitatively understood in physical terms, 
namely the idea that the evolution of the system before freeze-out is dominated by the longitudinal motion. Thus, 
the hydrodynamic transverse motion can be neglected or at least factorized out in order to study the longitudinal 
flow only. 

Indeed, the two seminal applications of relativistic hydrodynamics to particle and heavy-ion collisions, by Landau 
[3] and by Bjorken [4], start with this 1 + 1 approximation, valid in the determinative stage of the reaction. The 
longitudinal hydrodynamic approach has found many applications. It has been used in the literature [5, 6] in order to 
discuss aspects of the hydrodynamical flow which are relevant for the physical understanding of high-energy particle 
' scattering and, more recently [7], of heavy-ion collisions. 

Soon after the first proposal by Landau and its derivation of a large-time approximation [3] , Khalatnikov [8] showed 
that (1 + 1) hydrodynamics derive from a potential verifying a linear equation. The Khalatnikov potential has been 
used in the literature in an initial period [9, 10], but has not been recently considered, to our knowledge. Very recently, 
the interest on looking for exact solutions of (1+1) hydrodynamics has been revived and one finds new examples and 
applications of exact solutions, e.g. [11, 12]. For instance, in a recent paper [11], a unified description of Bjorken 
and Landau (1+1) flows has been proposed as a class of exact solutions of (1+1) hydrodynamics based on harmonic 
flows. (1+1) hydrodynamics appears also quite recently in the application of string-theoretical ideas to the formation 
a strongly interacting quark-gluon plasma [13]. These exact solutions allow to find explicit analytical solutions for 
relevant observables. Among these, the entropy flow dS/dy, where y is the rapidity, is quite interesting, since it may 
be related to the multiplicity distribution of particles. Our goal is to go beyond particular cases and obtain a general 
expression of the entropy flows as a function of the temperature for a generic solution of (1+1) hydrodynamics, i.e. 
for a generic solution of the Khalatnikov equation. 

We are interested in the distribution dS/dy of entropy density per unit of the rapidity y which is related to the 
flow velocities u± — e ±y in the (1+1) approximation. This "hydrodynamic observable" depends in an essential 
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way on the assumed hypersurface 1 through which we want to compute the flow (and eventually relate it to physical 
observables in collisions). For given entropy s and energy e densities, dS/dy depends on the hypersurface one considers 
to follow the hydrodynamic evolution. It is particularly interesting to consider hypersurfaces corresponding to a fixed- 
temperature. Indeed, it allows to follow the cooling of the hydrodynamic flow, from an initial stage characterized by 
a high temperature, towards a final stage which is often associated to a freeze-out temperature. Hence our aim is to 
derive an expression of dS/dy as a function of temperature for a given Khalatnikov potential and to investigate its 
properties, both on theoretical and phenomenological points of view. 

Our main new result, i.e. the general expressions for the entropy flow as a function of the Khalatnikov potential, 
can be found in three different versions (49, 50, 53). One or the other expressions can be more suitable for a given 
explicit problem. 

The plan of our paper is as follows. First, in section II, we group, for completion, all the necessary material, 
including the hydrodynamic equations, the Khalatnikov potential, its equation and solutions, recast and derived in 
a modern framework using light-cone variables. In section III, we formulate and derive the general expression of the 
entropy flow in rapidity as a function of the temperature evolution. In section IV, we derive and study a family of 
exact solutions, namely the ones where the final entropy distribution is dominated by the hydrodynamical evolution 
and not by the initial conditions. They generalize the Landau flow, and give the asymptotic behavior of physical flows 
in the limit of long hydrodynamical evolution. We provide in particular the exact analytic expression of the final 
entropy distribution corresponding to the Belenkij-Landau [9] solution. Then, in section V, we compare the profile of 
the entropy distributions, as well as their energy dependence, with the relevant experimental data. The final section 
is traditionally devoted to conclusions and outlook. 



II. (1+1) RELATIVISTIC HYDRODYNAMICS OF A PERFECT FLUID 

A. Hydrodynamic equations 

We consider a perfect fluid whose energy-momentum tensor is 

T"" = (e + p)u V - pif v (1) 

where e is the energy density, p is the pressure and (^ = {0, 1, 2, 3}) is the 4- velocity in the Minkowski metric rj^ u . 
It obeys the equation 

d^ v = . (2) 

Using the standard thermodynamical identities (where we have assumed for simplicity vanishing chemical potential): 

p + e = Ts ; de = Tds ; dp — sdT , (3) 

the system of hydrodynamic equations closes by relating energy density and pressure through the general equation of 
state 

We consider now the (1+1) approximation of the hydrodynamic flow, restricting it only to the longitudinal direction. 
Within such an approximation, the effect of the transverse dimensions is only reflected through the equation of state 
(4). Note that we do not a priori assume the traceless condition T w = 0, and thus the fluid is considered as "perfect" 
(null viscosity) but not necessarily "conformal" (null trace). 
Let us introduce light-cone coordinates 



z 



where r = Vz + z~ is the proper time and rj = | ln(z+/z ) is the space-time rapidity of the fluid. We also introduce 
for further use the light-cone components of the fluid velocity 



= u° ± u 1 = e ±y , (6) 



We keep the term hypersurface in the (1+1) case to keep track, even formally, of the transverse motion. 
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where y is the usual rapidity variable (in the energy-momentum space). 
The hydrodynamic equations (2) take the form 

(d + + d-)T m + (d + -d-)T 01 = 

(d+ + d-)T 01 + (d+-d-)T 11 = 0. (7) 
B. Khalatnikov potential 

It is known [8, 9] that one can replace to the non-linear problem of (1+1) hydrodynamic evolution with a linear 
equation for a suitably defined potential. In this section we follow the method of Ref. [8] , recasting the calculations in 
the light-cone variables. 

Inserting in (7) the known relations (1) for T^ u &nd expressing everything in light-cone coordinates using (5,6), one 
obtains the following two equations: 

e 2y - 1 \ / 1 - er 2y \ 

-^—)d + (e+p)+e 2y (e+p)d + y+ I J cL (e + p) + e~ 2y (e +p) d-y + 8+p - d-p = 

S d+{e+p)+e 2 y{e+p)d+y+ I — j d_ (e + p) - e~ 2y (e +p) d.y - 8+p - d.p = . (8) 

In (8) the energy density e and pressure p are considered as functions of the kinematic light-cone variables (z + , z~ ) . One 
key ingredient of the potential method [8] is to express the hydrodynamic equations in terms of the hydrodynamical 
variables y — logu + = — logu~ and 9 = log [T/T ] , where T is an arbitrary temperature scale. 

Relations (8) can be further transformed by inserting the differentials of the thermodynamic relations (3), namely 

d ± (e + p)=T Q d ± ( S e e ) 

d ± p = T s d ± e e . (9) 

Multiplying the first equation of (8) by (e~ 2y + 1), the second by (e~ 2y — 1), adding them and using (9), one obtains: 

8+ (e e+y ) - 8- (e e - y ) . (10) 

Eq.(10) proves the existence of a potential 2 $(z+,z _ ) verifying: 

d T ^(z + ,z-) = u ± T = T e e±y . (11) 

In this way, (10) is automatically satisfied. 

In order to transform the system of equations (8) from the kinematic variables (z + , z~) to the dynamical ones (9, y), 
one introduces [8] the Khalatnikov potential x, considered as a function of (m + T, u~T) through a Legendre transform: 

X(w + T, vTT) = $(z+, z~) - z-u+T - z + u~T , (12) 

where z^ are functions of (u + T,u~T) implicitly defined by (11). Hence, we get: 

j,*, = + [a + * - u-T] J*^ + [«.. - u*T\ ^ s , (13, 

where, due to the relations (11), the terms between brackets are zero. Knowing the Khalatnikov potential \i which 
is a function of the thermodynamic variables, one can find the kinematic variables of the flow by derivation. 

In the following, we will always consider the Khalatnikov potential x as function of 9 and y, keeping the same 
notation for \- That change of variables corresponds for the differential operators to 

e-°T y (d e ± d y ) . (14) 



d{u±T) 2T 



2 The function <I> has some degree of arbitrariness since we could define = Toe s ^ y + <p^(z^), with <fi—(z ) and (^_|_(z+) arbitrary 

one-variable functions. This freedom, analogous to a gauge choice, does not modify the final results. 
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In those variables, relation (13) writes 

zH0,y) = e~ e± y (-d eX ±d v x) ■ 
From (15), one also gets the expressions for the proper time t and the space-time rapidity r\ (defined as in (5)) 

r(e,y) = ^ yW) 2 - (Ax) 2 



(15) 



v(0,y) 



2T 

y + 2 lo ? 



-Oex-OyXJ \dexj 



(16) 



C. Khalatnikov equation 

Coming back to the system of equations (8), another independent combination can be obtained. Multiplying the 
first equation by (e~ 2y — 1), the second by {er 2y + 1) and adding, we obtain d + (e v s) + d- (e~ y s) — 0, or equivalently: 



d+ (u+s) +d- (u~s) = 0. 



(17) 



This relation corresponds physically to the conservation of the entropy along the flow. It is a property of the perfect 
fluid that the motion of the pieces of the fluid along the velocity lines is isentropic. 

Following the logics of the Legendre transform, we transform relation (17) using the (8, y)-base. For this sake, we 
write down the following partial derivatives: 



9(u±s) _ 



. ds 

de 



dy 



= ±u*s = 



9(u ± s) dz + 9(u ± s) dz 

dz+ df + dz~ df 
9(u ± s) dz + 9(u ± s) dz~ 



+ 



Solving this system of linear equations we obtain: 



d(u + s) 


1 


+ dz ds 


dz+ 


D 


U dy d9 


d(u~s) 


1 


_ dz + ds 
U dy dO 


dz~ 


~D 



dz+ dy dz dy 



, dz 



(18) 



+ U 8- 



dz+ 



where' 5 



D 



dz + dz dz + dz 



de dy dy d9 
Inserting (19) into the entropy-flow conservation relation (17) we acquire: 



ds 

de 



_ dz 
dy 



dz+ 



dy 



_dz 

~df 



dz+ 

de 



o. 



Obtaining the expression of the derivatives from (15), the equation (21) leads to: 



1 ds 



[dex ~ d 2 y x] ~ dex + dgX = • 



s de l " u/v y/ 

Making use of the sound velocity relation (4) we finally arrive at the Khalatnikov equation [8, 9]: 

c 2 d 2 eX {e, y) + [l- c 2 ] dg X (e, y) - d 2 yX (o, y) = o. 



(19) 



(20) 



(21) 



(22) 



(23) 



Hence, the non-linear system of equations which governs the (1+1) hydrodynamical flow has been converted into a lin- 
ear, second-order, hyperbolic partial differential equation. Note that the Khalatnikov equation is valid independently 
from the specific form of the sound velocity. 



3 We assume that the determinant D is non-zero, which is the case except on exceptional lines [9] . 
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D. Application: solutions of the Khalatnikov equation for fixed c, 



In this section, for our purpose, we present the solutions of the Khalatnikov equation with a constant speed of 
sound: 

= f = a' 1 , (24) 

where g will be considered as a parameter in the Khalatnikov equation (23). Note that in this case the general 
relations (3) write 

e = 9P=e (|r) 9+1 = e e^, (25) 



for the energy density and 



for the entropy density. 
Writing 



and inserting it into (23), we acquire: 



X(9,y) = e-( £ ^) e Z(6,y) (27) 



d%Z - g tfZ - Z = 0, (28) 



where we have used a compact notation for partial derivatives. 

It is convenient to replace the variables 9 and y by a and 0, defined by 

a = -9 + — and = -0 - — , (29) 

such that equation (28) takes the form 

d a d p Z(a, 0) - Z(a, 0) = . (30) 

We solve this equation following the Green's functions formalism, i.e. we look for distributions G(a, 0) such that 

d a d p G(a, 0) - ^ G(a, 0) = S(a)S(0) . (31) 
The relevant solution of equation (31) is 4 

G(a,0) = 0(a) 0(0) I (i-1^0) , (32) 



2 

with Jo the modified Bessel function of the first kind and 6 the Heaviside function. Using the relation 

5(a) 6(0) = s(-9+^j6 (-6 - = ^g 5(9) 5(y) , (33) 
we deduce from (32) the relevant Green's function of (28): 



G(e,y) = -^-G(a,0) = -^-e(-e + ^-) e(-e-^-) iJ g —±Je*- v -\, (34) 



4 There exist other Green's functions of equation (31), with e.g. Q(—a) instead of ©(a), or ©(— 0) instead of ©(/3). Assuming that 
the fluid naturally expands and cools down during the evolution, and taking the arbitrary temperature scale To to be the maximal 
temperature of the sources, —8 increases with time. Thus (32) gives the only physical solution of equation (31), analogous to the 
retarded propagator of the D'Alembert equation. Finally, note that for the obvious physical requirement of finite behavior at a, (3 — » oo, 
we reject the solutions of (31) containing the Bessel- Kq function instead of Iq. 
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which verifies 

d 2 e G - g d 2 y G ~ {^L f G = 5{9)8{y) . (35) 
Thus, we can construct the general solution of Khalatnikov equation (23), inserting a distribution of sources F(6,y), 



as: 



X (0,y) = e-(^) e J dy J d6G(9-6,y-y)F(6,y) 



'^"'«r SF^ m J^L-if-^\. (36, 



4 Vff J Je+\y-y\/^9 

Equation (36) gives the most general solution, for any distribution of sources of hydrodynamic flow. In the context 
of heavy-ion collisions, we are mostly interested in solutions that correspond to the evolution of a flow starting from 
initial conditions on a curve of the (0,y) plane. Therefore, one should impose constraints on F(6,y), in order to 
decribe the initial conditions. In section IV we will consider a physically-interesting sub-class of solutions. 

III. DERIVATION OF THE ENTROPY FLOW 

Coming back to the general formalism, let us now derive the exact formula for the entropy flow dS/dy at a given 
fixed temperature Tp = Tq e Sp , as a function of rapidity y. For a general (1+1) hydrodynamic expansion we consider 
the solution formulated in terms of the general Khalatnikov potential %(#, y), given by (36). The entropy distribution 
at fixed temperature is expressed through the amount of entropy flowing through the hypersurface of fixed temperature 
Tp, in an infinitesimal rapidity interval. It is given by (see e.g. [11]) 

dS u»d\„ „ dX 

— = sp—^- = s F u^— , 37 
dy dy dy 

where d\ is the infinitesimal (space- like) length element along the hypersurface of fixed temperature Tp, and n M is the 
normal to the hypersurface. The entropy density depends only on the temperature and not on y. Hence it is constant 
along the fixed-temperature hypersurface, namely sp = s(Tp) oc Tp . 



A. The flow through the fixed-temperature hypersurface 

As we have mentioned, we concentrate on hypcrsurfaces at fixed temperature Tp (or equivalently at 9p = 
log [Tp/To]). It is convenient to use as kinematical functions the proper time r = \J z + z~ and the space-time 
rapidity r\ = iln(z+/z~), considered as functions of 6 and y. In this (#,y)-base, the fixed-temperature hypersurface 
is parameterized by 

TF(y) = T(6 F ,y) 

vf(v) - v(0 F ,y) (38) 

considered as functions of y at Of fixed. The tangent vector to the hypersurface reads: 

V+(y) = z+'(y) = {r' F +r,'pT F )e^ 

V- (y) ee z F ' (y) = (r' F - n ' F T F ) , (39) 

where the primes denote derivatives with respect to y. Hence, we can construct the normalized perpendicular vector 
to the fixed-temperature curve (n + (y) , n~ (y)) defined by 

n + (y) n~(y) = 1 

\ [n+(y) V-(y)+n-(y) V+(y)} = . (40) 

Using (39) the second equation translates into 

n+(y) e-^ ( V ' f t f - t' f ) = n"(y) e" F {rf F r F + t' f ) . (41) 
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Provide \rf F {y)\ > for all y, we find 





r' F 


Vf t F ~ 




Vf t f - 




Vf t f + 





n + (y) = 
n~(y) = 

Following Ref.[ll], dX^ = dX is denned such that 

(dX) 2 = dX^dX^ = -dz+ dz F = -(r'p 2 - t f VF){dy) 2 , 
where the minus sign comes from the fact that the hypersurface is a space-like curve. Thus, we have 

dX = \Jrpip 1 - Tp dy . 

So, inserting (42) and (44) in (37), we finally find 

dS , 



dy 



(y) = s F [r F (y) rj' F (y) cosh(r] F (y) - y) + r' F (y) sinn^y) - y)} . 



(42) 



(43) 



(44) 



(45) 



B. Expression of the entropy flow 

Let us now introduce the expression of the entropy flow in terms of the Khalatnikov potential. Starting from (15), 
we obtain: 



cosh(?7 — y) = — 
sinh(?7 — y) 



1 



2rT e 
1 



9 dex{0, V) 

d v x{8, y) ■ 



2rT e e 

Now, inserting (46) in (45) we can eliminate the hyperbolic trigonometrical functions acquiring: 



dS x \-My) V'pjy) dex{0,y)+r'p(y) d yX (9,y)] 

dy [y> SF 2T e^T F (y) 



(46) 



(47) 



e=e F 



Furthermore, by differentiation of the relations (16) with respect to y, at 6 = 9 F , we find 

e- e [(dex)(d y dex)-(dyx)(d 2 y x)} 



My) = 
My) = 



2 Vidox) 2 - (dyX) 2 

[(dex? ~ (dyX) 2 + (dyX)(d y dox) - (dex)(d 2 x) 



(dex) 2 ~ (d v x) 2 



i=e F 



Then, inserting the relations (48) and (16) in (47), we obtain a remarkably simple expression, namely: 

dS . . sp 
dy™' = 27> 



^(y) = ^[d 2 y x^y)-9ex(o,y)]\ d=eF . 



(48) 



(49) 



which possesses a full generality, as long as the Khalatnikov potential x(# ; y) exists. In addition, using the Khalatnikov 
equation (23), (49) can be also written as: 



dS s F c 2 (T F ) 



dy 



(y) 



2T F 



[d 2 gX (e,y)-dex(0,y)}\ f)=B 



(50) 



There is an interesting third version of equations (49,50), featuring the potential $ instead of x- The definition 
(12) of x can be written alternatively 



X = $ - Tre n ~ v - Tre^+y = $ - 2Trcosh(?7 - y) . 



(51) 
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Inserting in (51) the first relation of (46), one obtains 

$ = x (0,y) - d e x(0,y) ■ (52) 

Inserting that last relation into (50), and considering the potential <& (originally defined in (11) as a function of z + and 
z~) now as a function of 9 and y, one gets a third equivalent formula for the entropy flow through fixed-temperature 
hypersurfaces, namely 

« W __2|S>) W ,,„ U . (53) 

The set of expressions (49), (50) and (53) form our main formal result. They provide the exact form of the entropy 
flow along fixed-temperature hypersurfaces, for a general (1+1) hydrodynamic evolution. We also mention that, 
beyond the derivation of the Khalatnikov potential at fixed sound velocity, formulae (49,50,53) still hold for a general 
speed of sound, once the solution of the general Khalatnikov equation (23) is known. It is important to note that 
relations (49,52) are valid as long as there exist x or $ potentials, even if there is no reduction to a linear equation, 
i.e. no entropy conservation in the (1+1) projection of the flow, while relations (50,53) are valid when the Khalatnikov 
equation holds i.e. with entropy conservation in the (1+1) projection of the flow. 



C. Examples 



Let us check the general formulae for the entropy flow considering exact hydrodynamical solutions known in the 
literature, namely the Bjorken flow [4] and the harmonic flows [11]. 

a. Bjorken flow: The Bjorken flow corresponds to boost-invariance, i.e. d y x = 0. In this case, the Khalatnikov 
equation (23) reduces to 

X "(e) + (g-l)x'(e)=0, (54) 

which has as generic solution 

X {6) = Ce-^ e , (55) 

C being an integration constant. Let us choose C and the arbitrary temperature scale To such that at the proper 
time r = r , the temperature of the fluid is T = T . Inserting (55) in relations (16) one finds C = 2T To/(g — 1), and 
the known expressions for the Bjorken flow, namely 

t(6, y) = To e~ 90 and r) = y , (56) 

i.e. the equality of rapidity with space-time rapidity. Finally, the Khalatnikov potential for the Bjorken solution 
writes 

Inserting (57) into (49), one obtains the entropy flow 

^r(y) = s f t f = sqt = est. . (58) 
dy 

Hence, as expected from boost invariance of the Bjorken flow, not only the total entropy but also the entropy flow is 
conserved. 

b. Harmonic flows: Following Ref.[ll], one is led to introduce new auxiliary variables l + {z + ) and l~(z~ ) satisfying 

£ - » ■ (») 
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where A = est. The thermodynamic variables can be explicitly written [11] as 5 

n 9 + 1 /,+2 , ,_2\ ,5-1 



49 r+n+'-w^ 



Using the property (11) of the potential $ one writes 

<9$ dz ± „ , cfe* 



Now, inserting (59) and the expressions (60), one obtains 

|| = AT e l± 'e e ™ = AT e W ++r f . (62) 
The expressions (62) are symmetric in l + and l~ and thus, by mere integration, one gets 

r l++r 

$(l + ,r) = AT / dv e~ v , (63) 
where the potential can be expressed in terms of 9 and y through 

i+ + r = V2\ y \ [- e -J 92 -j] ■ ( 64 ) 

Now, using our relation (53), and the relation 

de-f = AT e^ (l++r)2 d e (l + + P) , (65) 

one gets the result for the entropy flow 

| M =^^(«) (66) 

Using our general formalism, we thus recover the nontrivial result obtained by direct calculation (see Ref.[ll], formula 
(58)). Interestingly enough, we note that for the family of harmonic flows as an example, it appears to be much 
simpler to use formula (53) for the potential $ than using the Khalatnikov potential x itself. 

Note that we have to make a specific discussion of the limiting case when g = 1, that is when the speed of sound 
equals the speed of light. In fact in this case, the harmonic flow cannot be obtained as above and the solution for 
the flow acquires a more general form. Coming back to equation (23), one finds that the Khalatnikov potential itself 
is harmonic, namely x{^iV) = h+(y + 8^fg) + h-(y — 6^/g) where h+, h- are arbitrary functions. We thus recover the 
results noted in Refs.[15]. 



IV. EVOLUTION DOMINATED SOLUTIONS 



In general, a longitudinal flow in the final state follows from a longitudinal pressure gradient and/or from a longitu- 
dinal flow in the initial state. Let us consider the sub-class of solutions where the effect of the initial flow is negligible 
compared to the one of the initial pressure gradient. This sub-class corresponds to the dominance of the hydrodynamic 
evolution over the influence of the initial conditions. A typical example of such a solution is the Belenkij-Landau 
solution [9], where the fluid is initially at rest (the so-called "full stopping" initial conditions), and then expands into 
the vacuum. 



5 Here we use our convention 9 = logT/To with opposite sign w.r.t. [11]. 
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A. Khalatnikov potential and entropy flow 



In order to model an evolution-dominated flow, let us consider all the sources at rest, i.e. F(9,y) ax S(y). Let us 
also take the arbitrary temperature scale T to be the maximal temperature of the sources (hence 9 = log(T/T ) < 0), 
i.e. F(6,y) oc 6(— 6). All in all, we write 



F(9,y) = ^ K(9) B(-9) S(y) 
Inserting (67) in (36), and replacing the variable 9 by 9' = 9 — 9, one gets 

'.9-1 



x(0,y) 



f 

Je 



M 



v^' 2 - 2/7.9 K(9-9')d9', 



(67) 



(68) 



where the function K(9 — 9') carries the information on the initial conditions. Note that 9' is also negative. 

In the following it is convenient to use a Laplace representation of (68). Since 9 < 0, we introduce the Laplace 
transform, and its inverse, with respect to —9 as: 



Ki) = f d9 e< 6 f{9) 

J — oo 



/(*) 



70+400 d 1 _ lB 



5S ^ ™ • 



(69) 



where 70 is a real constant that exceeds the real part of all the singularities of the integrand, i.e the integral is 
calculated on an imaginary contour that lies on the right of all singularities. Following [9], the Khalatnikov potential 
(68) can be written as a convolution of the two functions : 



®{-e-\y\/Vg)h 



.9-1 



e(-9)K(9) = 



Vo 2 -y 2 /g 



70+ioo 

70 — ioc 
70+ioo 

70 — ioo 



2ni 

d"f 
2ri 



1 



i r~i 
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As the Laplace transform changes convolutions into ordinary products, one gets the Laplace representation 

f 70+ioo 



rto+100 1 

*V> y) = J £ 



K{j) 



Y 



(g-1) 2 



(70) 



(71) 



Notice that, while the expression of the solution (68) restricts the phase-space domain in the interval \y\ < ~^fg 9, 
equation (71) may allow for an analytic continuation of the solution of the Khalatnikov potential outside this region. 
However the outside region may be different (e.g. with x = 0, as in [9]). 

Let us now investigate the properties of the entropy flow given by the solutions (71) of the Khalatnikov equation. 
Inserting the Khalatnikov potential (71) into the expression of the entropy distribution (50), one is led to the following 
formula: 



- / 70+400 dj __ M7+V ) 



,70+^oc ^ 
J70 — ioo 27TI 



( 7 + 5 /2) 2 -l/4l K( 7 ) - 



(g-1) 2 



(72) 



In formula (72), one may distinguish the kernel 



Q(-y,v) 



(g-1) 2 



Y 



(g-1) 2 



(73) 



driving the dynamical hydrodynamic evolution as expressed on the entropy flow, and the coefficient function 



Cf(l)= ( 7 +.9/2) 2 -l/4 K(y) 



(74) 



which encodes the initial conditions of the entropy flow. 
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B. Total entropy 

Since we have a well-defined relation (72) for the entropy distribution, it is easy to perform the integration over y 
and obtain the total entropy flux through the hypersurface with fixed temperature T ~Tp. 
Formally, (72) leads to: 



Jo 



7 2 _ (g-1) 2 



(75) 



where we took into account the O (— 9 — \y\/y^g) function present in (70). Indeed, the hydrodynamical flow is limited 
in the region inside this domain, with possible contributions on the boundary 9 = -\y\/ yfg (Riemann waves, see, e.g. 

We know that, by construction, the flow is isentropic and thus the total entropy Stot is conserved. In fact it is 
possible to show that the dominant part of the total conserved entropy results from the kernel (73) more than other 
sources such as the coefficient function (74) or the boundary Riemann waves. Hence the hydrodynamical dynamics 
dominate. For this sake, let us release for simplicity the boudary limitations of the integral over y. One writes 



Indeed, the complex integral is obtained through the singularities of the integrand, which can be due either to the initial 
conditions (through singularities of if (7)) or to the hydrodynamical dynamics (through the pole at 7 = (g — l)/2), 
or both. If the singularities of if (7) are situated at the left (resp. right) of the pole, they will be subdominant (resp. 
dominant) in the total entropy. Assuming a dominance of the hydrodynamic flow, we get the final result of (76). The 
physical meaning of (76) becomes clear when using the thermodynamical relation (26) and the entropy density so at 
the temperature T . The total entropy writes 6 

S tot * S -^-K[(g-l)/2] , (77) 

and thus does not depend on the features of the flow at T = Tp. In conformity with the isentropic property of the 
flow, the total hydrodynamic entropy of the perfect fluid should be conserved, as the expression (77) is independent 
of Tp. This provides a self-consistency check for an evolution-dominated flow. In more general cases, one should also 
take into account the other contributions. 

A final comment is in order. A priori, the domain of integration \y\ < Y/2 comes from energy-momentum conser- 
vation. However, the formula (68) for the Khalatnikov potential is only valid in the domain \y\ < —^fg 9p. In the 
flow-dominated approximation, —^jOf and Y/2 are considered large enough such that the integration domain can 
extend to infinity and the kernel only contributes 7 . 



C. Relation to the Belenkij-Landau solution 

We have studied the dependence of our results on the coefficient function (74) by imposing various relevant analytic 
forms for if (7). We observed that typical meromorphic functions bounded by a constant 8 at infinity and with poles 
at the left of 7 = give smooth and similar entropy flow distributions, almost identical at large enough 9p. Hence 
we conjecture that all physical evolution-dominated solutions to be almost identical, at least for a sufficiently large 



Note that K is dimensionless as the potential x- 

We have performed numerical checks which show that thanks to the decreasing exponential behavior, the boundary term contributes 
negligibly to the total entropy. 

Indeed, choosing if (7) of strictly positive degree leads to an unphysical angular point at y = and to a function K(9), see (70), 
containing derivatives of Dirac distributions, i.e. structures that are too singular to describe physical flows. 



12 



value of Op — log(Tp/To). In order to provide an analytic expression for the entropy now characteristic of the family 
of solutions, we remark that the following choice of the coefficient function (74) 



C f {l) = C 7 + 



.9 



1 



k{ 1 ) = 



c 



(78) 



where C is a dimcnsionless constant, corresponds to the hydrodynamical flow with an initial full stopping condition 
[9]. 

The Belenkij-Landau solution [9] describes the evolution of a slice of fluid of width 2L initially at rest and expanding 
in the vacuum. It consists in a hydrodynamical flow bounded by Riemann waves. The matching conditions between 
the flow and the waves in space-time translated in terms of temperature and rapidity variables are realized by imposing 
zero boundary conditions on the Khalatnikov potential x on the characteristics = ±y/y/g. Another condition on 
the potential is that the center of the slice remains by symmetry at rest (y = 0) during the evolution. 

We have checked that the energy flow, resulting from modifications of the Ansatz (78) satisfying the dominance of 
the kernel singularity, is not sensibly modified from the one given by inserting the coefficient function (78) into (72). 

Inserting now (78) into (71), the Khalatnikov potential between the characteristics —0 > \y\/y/g acquires the 
analytic form [9, 10] 



x(0,y) = C 



•/a 



9- 



l0'2 



(79) 



The potential is identically zero in the region ~0 < \y\/y/g. Note that the constant in (78,79) is such that C oc LT 
with our notations. 

Let us now insert this specific solution to our general formula (50) for the entropy distribution. Calculating the 
derivatives we find: 



dS_ 

dy 



(y) = s F 



(9-i)c 

AgT F 



(a-i) 



Of 



9-1 



y 2 /g h 



.9-1 



4 



y 2 /g 



Of 



<M y 2 /9 



(80) 



Since Of is negative, this expression is always positive, at least in the region 2 F — y 2 /g > 0. Hence the positivity of 
the entropy flow is ensured. Finally, the expression (80) is divergence-free, since it is finite for 9 2 F ~ y 2 / g. We also 
note that it is still real in the analytic continuation of the solution (80) for 9 2 F < y 2 /g, since in this case both the 
numerator and the denominator of the second term are purely imaginary 9 . 

For the total entropy, inserting (78) into the general formula (75) one gets, using the thermodynamic relations (26) 



Stot — 



C s F 
V9T F 



-{ g -i)e F _ 



C sq 

V9T0 



(81) 



A comment is in order at this point. Condition (78) has been considered to describe the so-called full stopping 
conditions. In the original papers [9], it consists in the assumptions that i) there is a specific plane where the medium 
is at rest for all times, and that ii) on the vacuum-boundary we have just a simple (Riemann) wave. In fact, the 
resulting entropy flow distribution is expected to be more general and is characteristic of the evolution-dominated 
hydrodynamic flows. Hence the Khalatnikov potential (79) (already obtained in [9]) and the entropy flow (80) may 
serve as an analytic formulation for the class of evolution-dominated flows. In fact, their features are essentially 
determined by the evolution kernel Q(7, y) (see formula (73)). 

Finally, it is interesting to note that formula (72) gives the possibility to compare the hydrodynamic predictions 
with those of other existing models of heavy-ion (and eventually hadron-hadron, soft scattering) reactions. This 
relies on the possibility of relating thermodynamic quantities, such as the temperature and the entropy, to observed 
properties of the particle multiplicities. In this scheme, the rapidity of particles is defined by the corresponding 
value of y = logu + obtained from the fluid velocity of the lump of fluid giving rise locally to the hadrons. In the 
same context, the overall temperature gradient Of will be related to the total available rapidity and the entropy to 
the multiplicity up to phcnomcnological factors. We will discuss in the next section the phcnomenological issues of 
our derivation, but the general theoretical idea is that formula (72) can be compared with the one-particle inclusive 
hadronic cross section which is related to the scattering amplitudes. In this respect the generic form (73) of the 
hydrodynamic evolution kernel may serve for a comparison of hydrodynamic properties with conventional models of 
scattering amplitudes. 



Positivity may also extend, but is not ensured due to the appearance of Bcssel function zeroes. 
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V. PHENOMENOLOGICAL APPLICATIONS 



Motivated by the seminal works of Landau [3] and Bjorken [4], the comparison of their predictions for (1+1) hydro- 
dynamics with some features of the data has been made (see e.g. [5-7, 9, 10]). Even if such a rough approximation, 
ignoring the details of the transverse motion or of the hadronization, cannot replace the numerical simulations, it 
has given some useful information on the dynamics of the quark-gluon plasma. For instance, the order of magnitude 
estimates made using the Bjorken flow in the central region [4] and the comparison of the multiplicity distributions 
with the predictions of the Landau flow [6, 7] have indicated that the proper-time region during which the hydrody- 
namic flow is approximately (1+1) dimensional has a deep impact on the whole process. Our aim is to take benefit of 
the explicit form (80) representative of the entropy of evolution-dominated flows, based on the Khalatnikov potential 
(79), to revisit the discussion in the light of recent experimental results. 

For the phenomenological application, we will concentrate on the entropy flow corresponding to the Bclcnkij-Landau 
solution. From (80) and (81), one obtains the formula 



dS (g - 1) <£^ii ( 

-T-{y) = Stat . ' e ' 



.9-1 



y 2 /g 



.9-1 



y 2 /g 



y/o'i 2/7.9 



(82) 



where we have used the normalization by the total entropy Stot- Formula (82) still depends on two hydrodynamical 
parameters F , the logarithmic temperature evolution, and the speed of sound c s = g^ 1 ^ 2 . 



A. Multiplicity distribution at fixed energy 



In order to investigate the phenomenological validity of formula (82), let us consider the experimental BRAHMS 
data for the charged multiplicity distribution in the most central collisions as a function of the rapidity measured 
recently at RHIC [16]. For sake of simplicity, in accordance with the (1+1) dimensional approximation of the dynamics 
that we consider, we will make the following assumptions. We will identify the rapidity if the fluid elements yf = 
l/21og(w + /u~) with the rapidity of the particles y p = 1/2 log(p+/p~). We thus keep the same notation y. In the 
same way, we assume that the multiplicity distribution of produced particles dN /dy in rapidity can be considered to 
be equal, up to a constant factor, to the entropy distribution 10 dS /dy. One expects that the end of the hydrodynamic 
behavior appears at a typical temperature Tp, related to a hadronization or freeze-out temperature, and independent 
of the total c.o.m. energy of the collision. On the other hand, the initial temperature Tq is expected to depend on 
the total c.o.m. energy (or equivalently on the total rapidity Y) and on the centrality of the collision, through the 
energy density e(To) of the medium produced by the pre-hydrodynamic stage of the collision. Thus, Of — logTp/To, 
should be a function of Y and of the centrality. Our formalism, based on the 1 + 1 dimensional approximation of the 
flow, is not appropriate to have a precise description of the freeze-out. Note however that some improvement could 
be obtained by using, e.g., the Cooper-Frye formalism [17] in the derivation of the entropy flow. We postpone this to 
further studies. 

Using then formula (82) for dN /dy and fitting BRAHMS data by adjusting the parameter Of we obtain a good 
description for different values of g. In Fig. 1, as an example, we present the BRAHMS data fitted with (82), for four 
pairs of g and Of values, reported on the figure. In these plots the solid line corresponds to the physically meaningful 

2 2 

region (^- < Op), while the dotted line corresponds to the analytic continuation of formula (82) in the region ^- > 6 F , 
where the applicability of the solution (82) is theoretically questionable. 

The phenomenological application appears to be correct for quite different values of the speed of sound c s = g^ 1 ^ 2 . 
The overall form of the curves is satisfactory. For the first curve at c s = 1/V3 (i.e. the conformal case), however, the 

analytic continuation beyond y < 2 F is soon reached 11 . We will comment on this remark later on. Indeed, when 

decreasing the speed of sound, e.g. for g = 5, the physical domain ^- < 2 F extends in rapidity. 
Some comments on these results are in order. 

a) It has been well-known since long [6] and confirmed more recently that a Gaussian fit to the data 

f(y)~e-y 2 / Y (83) 



We also assume that the multiplicity distribution of charged particles is proportional to the total one. 
One may also note that the curve indicates a violation of positivity before the kinematical limit. 
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FIG. 1: BRAHMS data fitted with the hydrodynamic formula. The data are taken from [16] and they correspond to charged 
pions in cental Au+Au collisions at ^s NN = 200 GeV. The solid line corresponds to the physical region (^- < Bp), while the 

2 2 

dotted one corresponds to its analytic continuation (^- > F ). The small vertical line marks the experimental beam-rapidity 

Yheam ~ V/2. 



with a variance W, as predicted by Landau [3], was reasonably verified. We noticed, that, indeed, expression (82) 
has an approximate Gaussian form, but it does not correspond, except for very large 9p, to the expansion of the exact 
entropy distribution near y = 0, as in the original argument [3] which was based on an asymptotic approximation 12 . 
Hence, the sub-asymptotic features of the full solution plays an important phenomenological role. 

b) There is apparently no track of the transition between the physical regime ^- < 6 F and its analytic continuation, 
described by the dotted lines in Fig.l. This is related to the mathematical property of the general solution (72) 
expressed using a Laplace transform. In short, the 7 ,i Bessel functions are transformed into J ,i with the same 
argument up to a factor i, without discontinuity. 

c) This transition is however meaningful. In fact one knows that the lines = ^f6 F delineate different regions of 

the hydrodynamical regime. Discontinuities, and thus shock or Ricmann waves may occur at these boundaries, called 
characteristics of the equation [18]. Hence some other solutions may branch at this point (see e.g. [9, 10]). However, 
our results do not depend on the specific form of these other solutions such as the shock waves considered in [9, 10]. 



B. Energy dependence of the multiplicity distributions 

Going a step further, we would like to interpret the energy dependence (i.e. the Y dependence) of the (1+1) solution 
for the entropy flow compared with multiplicity data. For this sake, we make use of the Gaussian fits reported 13 in 
[16] for different sets of data ranging from the AGS to RHIC. 

In Fig. 2 we give the determination of ^ as a function of Y which gives a good description of the Gaussian fits 
with the variance taken from [16]. As in the previous study of BRAHMS data, we performed this fit for six different 
values of g. As shown in Fig. 2 the corresponding relation is clearly linear. We can write: 

Of = -n{Y-Y ) . (84) 

As shown on the figure, the constant term Yq depends appreciably on g while the slope k 0.2 remains only 
slightly dependent on it. On a physical ground, the linear relation (84) has a reasonable interpretation. The initial 



It is indeed easy to verify that for phenomenological values of Bp, this approximation does not work in the data range. 

In fact, the prediction (83) fits reasonably well, but we used instead the actual best-fit determination of the variances provided in [16]. 
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FIG. 2: The hydrodynamical parameter 6p as a function of Y. We describe the dependence for for six values of g, in a large 
range starting from the canonical value 3 which would correspond to a conformal fluid in (3 + 1) dimensions. 



temperature of the medium is expected to grow as a power k < 1 of the incident energy. One finds approximately 
T /T F = e~ 6p ~ e -22(r-Yo), Hence the more energy is available, the longer the hydrodynamical evolution lasts. At 
smaller speed of sound the hydrodynamical evolution has to occur on a larger temperature interval in order to describe 
the same entropy distribution, as could be expected. 

Moreover, there is a physical argument, analogous to the one proposed by Landau [3], for the existence of a linear 
relation (84) between the temperature ratio and the total c.o.m. energy. Assuming the approximate validity of the 
Bjorken relation 14 tq/tf m (Tf/Tq) 9 , where To (resp. Tp) are the initial (resp. final) proper-times of the (1+1) 
hydrodynamical evolution and reporting in (84) one finds 

\ogT /T F ^ng(Y-Y Q ) . (85) 

Indeed, following [3], the separation proper-time from (1+1) hydrodynamics to the (1+3) regime is of order log At, ~ 
iy where A is the typical transverse size of the initial particles. Assuming that we can approximate At, by tq/tf, 
and taking into account the Bjorken flow approximation, formula (85) is suggestive of the proportionality property. 
We leave the precise values for k, t /t f and Yo to a further determination of g since the data we discussed do not 
prefer a precise value of ng (to be compared with 1/2 obtained for log Ar s ). 

An interesting consequence of the linear relation (84) between 9f and Y at fixed g is the possibility of relating the 
general hydrodynamic entropy distribution (72) to the one-particle inclusive cross-section and thus to the appropriate 
scattering amplitudes. These are not easy to formulate in the hydrodynamical formalism. Being more specific, let us 
transform the formula (72) in terms of the energy dependence using k as the coefficient of proportionality in (84). We 
get 



f(v,Y) oc r i °° *L e-.Cr-*)^) [ (7 + g /2) 2 -l/4] *( 7 ) °^f^l . (86) 

<ly Jjo-ioC 27TJ L J / 2 _ (g 1)2 

Formula (86) shows that the characteristic hydrodynamic kernel Q (see (73)) appears also as the kernel of the Laplace 
transform in Y of the one-particle inclusive cross-section, up to a redefinition of the conjugate moment u> = wry of 
the total rapidity Y. This relation may be useful to compare various theories and models for scattering amplitudes of 
high-energy collisions with the predictions of hydrodynamic evolution. 



This relation, properly stating, is exact only for the Bjorken boost-invariant flow. However, one expects that it remains approximately 
valid in the central region of more general flows (see, e.g. [11].) 
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VI. CONCLUSIONS AND OUTLOOK 

Let us summarize the results of our study: 

On the theoretical ground, we have the following results: 

i) We have recalled and reformulated the derivation of the Khalatnikov potential and equation in terms of light-cone 
variables. This allows to formulate the initial nonlinear problem of (1+1) hydrodynamics of a perfect fluid in terms 
of solutions of a linear equation. As an application, using the Green function formalism we derive the general form of 
the solution for constant speed of sound. 

ii) Expressing the flow of entropy through fixed-temperature hypersurfaces, we provide general and simple expres- 
sions of the entropy flow dS/dy in terms of the Khalatnikov potential 15 . 

iii) We check and illustrate the simplicity of the obtained formulae for dS/ dy by applying the formalism to some 
exact hydrodynamic solutions which were not using the Khalatnikov formulation, such as the Bjorken flow and the 
less straightforward example of the harmonic flows of Ref.[ll]. 

iv) We use our formalism to find the entropy flow for the subclass of solutions for which the hydrodynamic evolution 
dominates over the infliucncc of initial conditions. A characteristic example of such flows is the one studied long ago 
by Landau and Belenkij [9] , corresponding to full stopping initial conditions. We provide an exact expression for the 
related entropy flow. 

As a phcnomenological application, we discuss the relevance of the full stopping entropy flow for modern heavy-ion 
experiments which was advocated, e.g. in Refs.[6, 7]. 

i) The exact expression of dS/dy for the Belenkij-Landau solution, depending only on the ratio Tf/Tq and on the 
speed of sound c s , is in agreement with the shape of the multiplicity distribution of particles dN /dy(y,Y) observed 
in heavy-ion reactions, with a linear relation between the temperature ratio and the total rapidity log T n /Tp = 
k[Y-Y q (c s )}. 

ii) However, comparing our exact results with the asymptotic Gaussian predictions [3, 6] for the multiplicity 
distributions, we find that nonasymptotic contributions play an important role in the phenomenological description. 

iii) The speed of sound, which is the remaining parameter in our study, is not determined by the multiplicity 
distribution, since the phenomenological description seems satisfactory for a rather large range of the parameter 
g = 1/c 2 .. However, even if one does not see any sizable effect on the curve for dS/dy(y,T) one notices that the 
physical domain of the hydrodynamic expansion is restricted by the condition y 2 < g 9 F , especially for a speed of 
sound as large as the conformal one c s = 1/V3. 

This summary of conclusions leads to a few comments on possible further developments of our approach. Some of 
them are technical but could provide a further insight on the features of (1+1) hydrodynamics. First, it should be 
useful to study in detail a larger set of solutions. Second, implementing the Cooper-Frye formalism [17] directly in 
terms of the Khalatnikov potential could refine the hypothesis of a fixed final temperature Tp. Also, the investigation 
of the entropy flow through other hypersurfaces, in particular the proper-time ones (c/. [11]) would be welcome, in 
particular to allow for a straightforward implementation of fixed proper-time initial conditions. 

On the phenomenological point-of-view, it is important to develop the comparison of the (1+1) approach with 
the data and including more corrections to the idealized dominance of the longitudinal motion. One question could 
be settled at least phenomenologically, which is the determination of the best fit for the speed of sound parameter, 
which is presently rather free. Also, including a viscosity contribution is another important issue, together with the 
investigation of the entropy flows with varying speed-of-sound. 

One may ask what is the meaning of the transition on the lines y 1 = g9 2 Fl which appear even in the physical 
rapidity region. Mathematically, they are the Ricmann characteristics of the Khalatnikov equations and as such, they 
are regions where discontinuities may appear [18]. Indeed, these characteristics were used in the old studies [9, 10] 
to connect boundary Riemann waves to the domain of dynamical hydrodynamic evolution. What is their meaning, if 
any, in the modern view we have now of high-energy collisions, is an interesting open question. 

On a more conceptual point-of-view, our study of the entropy flow and its dependence on rapidity may have 
some impact on recent studies [13, 14] of the AdS/CFT correspondence. It relates the hydrodynamics of a fluid, 
whose microscopic description is the one of a gauge field theory, with the string theory in a higher-dimensional space 
where the Einstein equations govern the gravitational properties of its low-energy regime. The actual realizations 
of the duality correspondence for a collective flow require boost-invariancc and thus are limited to the Bjorken 
flow. This flow contains an infinite energy, and is thus of limited relevance. Knowing the analytic form of more 
physical solutions should be helpful to derive their dual gravitational backgrounds, which is a priori a formidable task. 



After the completion of this paper, we noticed a related study in Ref. [19] 
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